Vacuum polarization of a scalar field in wormhole spacetimes 
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I. INTRODUCTION 

o"- 

Wormholes are topological handles in spacetime linking widely separated regions of a single universe, or "bridges" 

joining two different spacetimes |jj. Interest in these configurations dates back at least as far as 1916 g] with 

punctuated revivals of activity following both the classic work of Einstein and Rosen in 1935 |3 and the later series of 

Q,| works initiated by Wheeler in 1955 Q. More recently, a fresh interest in the topic has been rekindled by the works of 

Morris and Thorne and of Morris, Thorne and Yurtsever 0. These authors constructed and investigated a class 

of objects they referred to as "traversable wormholes." Their work led to a flurry of activity in wormhole physics [M. 

One of the most central feature of wormhole physics is the fact that traversable wormholes are accompanied by 

t— I ' unavoidable violations of the null energy condition (NEC), i.e., the matter threading the wormhole's throat has to be 

possessed of "exotic" properties. The known classical matter does satisfy the usual energy conditions, hence wormholes, 

if they exist, should arise as solutions of general relativity and "unusal" and/or quantum matter .[] Recently, Barcelo 

and Visser flJIJlJ] have demonstrated the possibility of existence of traversable wormholes supported by classical non- 

minimally coupled scalar fields. In their investigation, scalar fields play a role of the "unusual" matter violating the 

NEC, and therefore, instead of the problem of the wormhole's existence, the problem of the existence of appropriate 

Q\ . scalars comes to the forefront. 

An alternative approach is to regard wormholes as semiclassical in nature, and seek them as solutions of the Einstein 
equations, using the expectation value of the stress-energy tensor as the source of gravity (self-consistent wormholes). 
Such an idea first appeared and was realized in Rcf. ||l5fl . Note that the main difficulty of this approach is that the 
expectation values have the strong functional dependence on the metric tensor g M „ and are generally impossible to 
O" 1 be calculated analitically. For this reason, one is forced to use approximate expressions for the expectation value 
5^ \ of the stress-energy tensor. In Ref. |15), the Frolov-Zel'nikov approximation pa] for (!),„) for conformally invariant 
bi). massless fields in static spacetimes, based on pure geometrical arguments and common properties of the stress-energy 
tensor, has been used. Later, in Ref. |17[ self-consistent spherically symmetric wormhole solutions have been found 
numerically in the framework of the analytical approximation for (T^„), which was derived by Anderson, Hiscock 
and Samuel pq| o n the basis of the WKB approximation for modes of a scalar field. Recently Taylor, Hiscock, 
and Anderson ]19|], using the Dewitt-Schwinger approximation, have calculated analytic expressions for the stress- 
energy of a quantized massive scalar field in a number of static spherically symmetric wormhole spacetimes. As well 
some arguments in favour of the possibility of existence of self-consistent wormhole solutions to semiclassical Einstein 
equations have been given in Ref. |2Cf| , 

Recently one of the authors [Ell has performed approximate calculations of (0 2 ) for a massive scalar field in static 
spherically symmetric spacetimes, which improve the result by Anderson, Hiscock, and Samuel |18|| . The improvement 
has been achieved due to the use of exact expressions for mode sums and integrals. Moreover, it has been demonstrated 
in Ref. [£1| that the result of approximate calculations depends on the appropriate choice of the zeroth-order solution 
of the mode equation. In this paper we apply the experience obtained in the previous work to derive an approximate 
expression for (4> 2 ) for a scalar field in a wormhole spacetime. 
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x The NEC can also be violated and wormholes can exist in theories based on modifications of general relativity, such as 
Brans-Dicke theory pHLOl, higher derivative gravity |llj or Gauss-Bonnet theory |12j. 



II. WORMHOLE SPACETIMES 

We shall regard a wormhole as a time independent, nonrotating, and spherically symmetric bridge connecting two 
asymptotically flat regions. The metric of wormhole spacetime (continued analitically into Euclidean space) can be 
taken in the form, suggested by Morris and Thornc 



ds 2 



f(p)dr 2 + dp 2 + r 2 (p)(d0 2 + sin 2 dip 2 ). 



(1) 



Here t = it is the Euclidean time, p is the proper radial distance, p € (—00, +00). We assume that the redshift 
function f(p) is everywhere finite (no event horizons); the shape function r(p) has the global minimum at p = 0, so 
that ro = min{r(p)} = r(0) is the radius of the wormhole throat. The topology of the wormhole spacetime is R 2 x S 2 . 
In order for the spacetime geometry to tend to an appropriate asymptotically flat limit at p — > ±00 we impose 



lim {r(p)/\p\} = 1, and lim f(p) = f±. 

p— »±oo p— >±oc 



(2) 



For simplicity we also assume symmetry under interchange of the two asymptotically flat regions, p <-> —p, that is, 

r (p) = r (-p) and f(p) = f(~p)- 

Note that introducing a new radial variable r instead of p by substitution 



r = r(p) 
allows to rewrite the metric (|lj) in the following form 

ds 2 = f(r)dr 2 + h(r)dr 2 + r 2 (d9 2 + sin 2 dip 2 ), 
with 



(3) 

(4) 



Hr) = -j, 



where the prime means the derivative with respect to p. 



III. APPROXIMATE EXPRESSION FOR 



2 ) 



Consider a quantized scalar field <f>; the scalar field is assumed to be both massive with the mass m and massless, 
and with an arbitrary coupling £ to the scalar curvature. We also assume that the field is in a vacuum state defined 
with respect to the Killing vector which always exists in a static spacetime. An unrenormalized expression for (<j> ) 
can be computed from the Euclidean Green's function Ge(x,x). Such the expression given in Ref. JL8) for the case of 
static spherically symmetric spacetimes reads 



(</> 2 (a;)) U nrcn = G E {x,T]X,f) 



1 r °° — 

= 732/ dwcos[w(T-f)]53 

J0 1=0 



(21 + l)p u iq u l - — yyj 



where the modes p^i and q^i obey the homogeneous radial mode equation 



d 2 S 
dp 2 



1 df 2dr 
2f dp r dp 



dS_ 
dp 



LU 2 1(1 + 1) 

f r 2 



The modes p u i and q^i also satisfy the Wronskian condition 



a 



dq u i dp^i 

PljI—, qwi 



dp 



dp 



r 2 Jl/2 ' 



m 2 + £R 



5 = 0. 



(5) 



(6) 



(7) 



2 The units % — c = G — 1 are used throughout the paper. 



where C u i is a normalization constant. 

Let us stress that points in Eq. ( [L9[ ) are splitted. Namely, points are separated in time so that e = (r — t), p = p, 
6 = 9, tp = tp. As was first pointed out by Candelas and Howard |2J| for the case of Schwarzschild spacetime, the 
Euclidean Green's function have superficial divergences with this separation of points. As discussed in Ref. Eq], these 
cannot be real divergences because the Green's function must be finite when the points are separated; to remove the 
divergences one has to subtract some additional counterterms. The terms (r/ 1 / 2 )^ 1 in brackets in Eq. ( |l9| ) are such 
counterterms (for details, see discussion in Ref. p8| ). 

Now let us proceed to the WKB representation by making the change of variables 



1 



Pul 



qui 



{2r 2 W) 1 / 2 

1 

(2r 2 W/) 1/2 



exp 



exp 



Wf-^dp 



Wf- 1/2 dp 



(8) 



where W is a new function of p. Substitution of Eq. (g) into Eq. (Q) shows that the Wronskian condition is obeyed 
if Cui — 1, and substitution of Eq. (g) into Eq. (0) gives the following equation for W: 



W 2 



I 



r'f r"f f'W fW" 3fW n 



^1(1 + 1) +mV + €/*+ ^ + -^ + tw + 2 liT 



2r 



2r 



4 W 2 



with 



R 



f" f' 2 2r'f 4r" 2/ ; 



/ 2/ 2 



rf 



(9) 



(10) 



Assuming that the functions / and r are varying sufficiently slowly, one can solve Eq. (S) iteratively with the small 
parameter 



^wkb = L+/L <C 1, 
where L is a characteristic scale of varying of the metric functionsn and 



u 



2 1 



-1/2 



(11) 



(12) 



Neglecting terms with derivatives in Eq. 



where 



we choose the zeroth-order solution as follows: 
W = Q, 



n 



I 



2C 



^1(1 + 1) + fm 2 + -| 



1/2 






./' 



r 



r 



Here we denote 



/i 2 =m 2 r 2 + 2(5 
with S = £ — 1/8. Stress that below we assume p 2 > 0. The fourth-order solution is 

W = n + T^ (2) + iy (4) 



(13) 



(14) 



(15) 



3 More exactly, one may define L as follows: L 1 = max -I [ln(/r)]' , [ln(/r)]"j , [ln(/r)]'"j , . . . >. 



with 



2 fi 4 ft 2 8 ft 2 8 ft 3 



(16) 



and 



W& 



W^ 2 



2<> 



/W^ /' y^ 2 )' f W'W^ f n'ww 3f n /2 w^ 
In 2 +_ 8~ n 2 ~~i n 3 8 n 3 + T ft 3 - 



(17) 



where 



V 



2r 



2r 



■Sf[R 



(18) 



Substituting the solution ( Jig ) into dq) and ( J19J ) , and neglecting terms of the sixth order and higher we can obtain 
the following expression for the fourth-order WKB approximation for (4> 2 ): 



Air 2 
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n 5 



i f 00 
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"16/ 
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(19) 



where 



t/=4M + 7T 



, , 2(5 
Mm 2 + ^ 



(20) 



A renormalized expression for (</> 2 ) can be obtained as follows: 

(0 2 )ren = Hm ((0 2 ) unre n ~ 



;ds 



(21) 



where the renormalization counterterms for a massive scalar field are given by 



;ds 



Gt>s(x,x') 

1 1 

8^V + 8^ 
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R 



C- 



ln 



'DSl 



16tH 



96tH 



r-R, 



a a aP 



a/3- 



(22) 



Here a is equal to one half the square of the distance between the points x and x' along the shortest geodesic connecting 
them. C is Euler's constant, R a p is the Ricci tensor and a a = ct ;q . The constant ttids is equal to the mass m of the 
field for a massive scalar field. For a massless scalar field it is an arbitrary parameter. A particular choice of the value 



of, 



'DS 



corresponds to a finite renormalization of the coefficients of terms in the gravitational Lagrangian. 



To perform the procedure of renormalization in practice we make use of the formulas obtained in Ref. |2] 
formulas give the expressions for mode sums and integrals as exact expansions in powers of e: 



These 



So(e,n) 



du cos 



,/i/a 



= 7^ + 2^ 2 -12»( ln 



1=0 



mf 1 ' 2 



2r 



1 + 1/2 



v/m 2 + n 2 + (I + 1/2) 2 



(") -^-^So(2nfi) + O(e 2 \ne), 



(23) 



ef 1 / 2 — - 
S n (e,n) = J du cos ( u ) \, 



ii 



a + 1/2) 2 



/■ / ^\u^ + ^ 2 + {i + i/2) 2 ] n + 1 / 21 



(2n-l)! 
n= 1,2,3,..., 

where (2n - 1)!! = 1 • 3 • 5 • • • (2n - 1) and 



2m - 1 ("- 1 ) ! ^ln^ + C7 > i + 5^, + 0( e 2 ln e ), 



2r 



(2n-l)H 



S'o(27r/i) 



x In 1 1 — x 2 1 



-da;, 



^(2^)=/ ^H--%U^ 



d /l d\ n ~ l s 2 ^- 1 ) 



(24) 



^Itt^ix i 1 * 



Together with the integrosums (|23[) , which diverge in the limit e — ► 0, the expression (19) contains the integrosums of 
the following form: 



N™(e,v) 



, ( tf 1/2 \^ q + i/gr 



n = 2,3,4,..., 

m = 1, 3, . . . , In — 3 



(25) 



The functions N™(e, /z) are regular in the limit e — > 0, and so one may directly set e = in Eq. (|25|). Then, changing 
the order of summation and integration in Eq. (|25|) and integrating over it gives 



^(Q,i*)= [^_i||| *TO, 



(26) 



with 



00 r; _i_ I V n 



,^> 2 + G+£) 2 ] 



Substituting Eqs. ( |23]j2q ) into Eq. ( |19| ) we may find the asymptotical expansion for (</> 2 ) U nren in the limit e — ► 0. 
Finally, earring out the procedure of renormalization (pl|), i.e. subtracting (</> 2 )dS; we obtain the renormalized 
expression for (0 2 ) in the framework of the fourth-order WKB approximation: 
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DS' 
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96/ 2 



-^(^ui 



i/V/ 



12/ I \r 2 ^ J 2 f Vr 2 ^ 
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IV. ANALYSIS 



The analytical approximation for (0 2 ), given by Eq. (p7|), may be separately analysed for two different cases: (i) 
r(p) > L, and (ii) r(p) -C £. 

A. The case r>L 

First, let us consider (c/> 2 ) in the region r(p)>L. Note that in this case /i 2 sw m 2 r 2 ^> 1. Really, it immediately 
follows from the basic WKB condition Xwkb = t*/£ € 1: 



assuming that £ is of order unity or less, we conclude that m 2 r 2 3> 1. The functions 5 n (27Tju) and N™{fx) have a 



simple asymptotical form at large values of mr; the respective formulas are given in Appendix (see Eqs. ( |A2| ), (A7)). 
Substituting these formulas into Eq. ( |27| ) and taking into account the asymptotical expansion ln(l + x) ~ x — ^x 2 
gives 

,2 2/ j2 n 5£-l/"" 5$ -Iff ( 1, 2 7, 13 \ /"/' 2 
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(28) 

The expression for {4> 2 ), given by Eq. (pi), has same order as the term proportionate to mT 2 in DeWitt-Schwinger 
expansion for (0 2 ), however it does not coincide with them. The reason of this is that the topological structures of 
Minkowski and wormhole spacetimes are different. 

As a simple example, we consider the wormhole spacetime with the metric 

ds 2 = -dt 2 + dp 2 + (p 2 + rl){d6 2 + sin 2 9 dip 2 ), (29) 

given by Morris and Thorne ||. Here f(p) = 1 and r{p) = (p 2 + r^) 1 / 2 . The characteristic scale of varying of the 
metric functions can be defined as follows: L~ 2 ~ r" Jr — ^o(/° 2 + r o)~ 2 near the throat, and L _1 ~ r' ' /r = p{p 2 +r 2 l )~ 1 
far from the throat. Stress that in both cases r(p) ~ L. Note also that the basic WKB condition L/L+ ^> 1 is satisfied 
provided r§jr c 3> 1, where r c — m^ 1 is the Compton length. Carrying out appropriate calculations in Eq. ( J28j ) we 
find 

4 7 r 2 m 2 (0 2 } rcn = *$> (30) 

where 

y M = - _ J na fix 4 - — x 2 - -£x 2 - -e + —) (31) 

u ; (x 2 + l) 4 V70 1120 4^ 2^ 21/ V ' 



and x = p/ro is a dimensionless proper radial distance. The obtained formula obviously shows that the qualitative 
behavior of (4> 2 ) is determined by the function ^ (x) and depends only on £. This dependance is very interesting. 
In Figs. 1,2 we present the family of curves of (cf) 2 } for various values of £. It is seen that (cf) 2 ) is purely negative at 
£ = 0, and becomes alternating for greater values of |£|. To characterize this feature more accurate we consider an 
averaged value of (</> 2 ), i.e., 

1 /"°° 1 / 8 1 0^7 \ 

_f) «*> - J^ L *« W V^^* ^ 3^S? fr * + 3« - iJSBo) • (32) 

It is easily to find that the averaged value is negative provided £i < £ < £2, where £1 « —0.7994 and £ 2 ~ 0.1744, and 
positive in the other cases. At last, we note that the average value of (cf) 2 ), computed in the metric (f29|), is negative 
for £ = (minimal coupling) and £ = 1/6 (conformal coupling). 

B. The case r < L 

Now let us analyse the analytical approximation (E7J) in the region r(p) <C L. Remember that L is the characteristic 
scale of varying of the metric, and so L" 1 is proportional to derivatives of the metric functions. Hence it follows from 



the relation r(p) <C L that one may neglect terms with derivatives in Eq. Ifflu with respect to terms without ones. 
As the result, we obtain 



4n A 



r 2 V 6 



In 



p 



T)S' 



-^ - ^ S {2tth) 



(33) 



It is worth to emphasize that an extreme example of the metric, for which the relation r(p) <C L is satisfied, is that 
when r{p) = r® = constant and f(p) = /o = constant, i.e., 



ds 2 



-fodt 2 + dp 2 + rl(d0 2 + sin 2 9 dip 2 ). 



(34) 



To illustrate a case of the mertic with non-constant coefficients obeying the condition r(p) <Lwe may consider the 
following example: f{p) = /o = constant and r 2 (p) = 7q[1 + a 2 ln(l + p 2 /po)], where a and po are some parameters; 
the relation r(p) <C L is satisfied provided ar /p a <C 1. 

Let, in addition to r(p) <C L, the relation r/r c 3> 1 be satisfied. Stress that it assumes that m ^ and so tods = to. 
In this case p « mr and we can simplify Eq. ( |33| ) by using the asymptotical form (Al) for Sq(2ttp) and taking into 
account that ln(/x 2 / (mDsr) 2 ) = ln(l + 26/(mr) T ) ~ 2S/(mr) 2 — 2S 2 / (mr) 4 : 



A7T 2 m 2 U 2 ) Icn 



— e- i +- 

2r 4 V 3 30 



(35) 



It is easily to see that the last expression is positive for all values £. 

Now consider the opposite relation mr = r/r c <C 1. Note that it permits in particular the massless case m = 0. In 
this case p rs {28) 1 ^ 2 the expression (p3) for (</> 2 ) can approximately be rewritten as 



4»vw = ^(e-J)'» 



2d" 



ml s r 2 



S 
2^2 



1 + 45 2nV25 



In particular, for conformal coupling £ = 1/6 (S = 1/24) we obtain <So(3 1 I 2 -k) w —0.0596 and 

0.0159 



4^ 2 (^ 2 ) re „ 



(36) 



(37) 



Finally stress that the obtained expression ( p3|) for (</> 2 ) is not approximate but exoci in the R 2 x 5 2 spacetime 
with the metric (p4J). 
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APPENDIX: ASYMPTOTICS FOR S N (2irp) AND N$f (p) 



In this appendix we obtain asymptotics at large values of argument for the functions S n (2np ) and N™(p). Consider 
the functions S n {2i:p) defined by Eqs. (Al). Making the substitution y = 2npx in Eqs. (|Al[) we find 



So(2irp) = A 5 



S n {2irp) 



dy, 



y\n\l — X 2 y 2 
ev + 1 



d ^-^Ty{yT, 



I (/ ^ n-1 y2(n _i) 

ev + l ! 



(Al) 



where A = (27r/i) _1 and /i = (m 2 r 2 + 2<5) 1 / 2 . Note that the integrands in Eqs. (Al) contain exponential functions 
and are exponentially decreasing at large values of y. Hence the main contribution into the integrals is provided with 
values of integrands in the region < y < 1. We are interesting in the case p^S> 1. In this case A<1 and Ay ^C 1 if 
< y < 1. Now we may use the asymptotical formula ln(l — A 2 y 2 ) = —X 2 y 2 — ^A 4 j/ 4 — ^A 6 j/ 6 + 0(A 8 y 8 ). Substituting 
this into Eqs. (Al) and integrating gives 
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1920(mr) 4 
S 2 (2xn) = 0((mr)- 4 ), S 3 (2x^) = 0((mr)- 6 ). 
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(A2) 



Consider now the functions 7V™(//): 



where n = 2, 3, 4, . . . and m = 1, 3, . 



and in addition. 
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2n — 3. It is easily to check that they obey the following recurrent formula 
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Using these formulas, one can express any function iV™(/i) via N^ifJ.). The function N^fi) can be represented as 
follows [see Ref. Pi, Eq. (5.1.26.23)]: 



iV 2 1 (u) = 

2 vw 4/x 



^' ' 2 _ * M ) ~^' ( 2 + * M 
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where V'( z ) — F'(z)/T(z) is the digamma function. Using the asymptotical properties of the digamma function (see, 
for example, Ref. p§|) and taking into account Eqs. (A4), (A5) we can obtain the following asymptotic for large 
values of //: 

N *M = 17T-T2 - 7^ + 0((mr)- 6 ), Nfo) = — I- - L^. + 0((mr)- 8 ), 
2(mr) 2 (mr) 4 4(mr) 4 (mr)° 
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FIG. 1. The curves in this figure display the the WKB approximation for (4> 2 ) (]28|) calculated in the large m limit in the 
Morris-Thorne metric (pa). From bottom to top the curves correspond to the values £ = 0, 1/16, 1/6, 3/10. 



FIG. 2. The curves in this figure display 
correspond to the values £ = 0, —0.1, —0.3, —0.6, —0.9 




} calculated for negative values of £. The labels a, b, c, d, e mark curves which 
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